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ABSTRACT 

We study the spectral evolution of PWNe taking into account the energy 
injected when they are young. We model the evolution of the magnetic field inside 
a uniformly expanding PWN. Considering time dependent injection from the 
pulsar and coolings by radiative and adiabatic losses, we solve the evolution of the 
particle distribution function. The model is calibrated by fitting the calculated 
spectrum to the observations of the Crab Nebula at an age of a thousand years. 
The spectral evolution of the Crab Nebula in our model shows that the flux ratio 
of TeV 7-rays to X-rays increases with time, which implies that old PWNe are 
faint in X-rays, but not in TeV 7-rays. The increase of this ratio is because 
the magnetic field decreases with time and is not because the X-ray emitting 
particles are cooled more rapidly than the TeV 7-ray emitting particles. Our 
spectral evolution model matches the observed rate of the radio flux decrease of 
the Crab Nebula. This result imphes that our magnetic field evolution model is 
close to the reality. Finally, from the viewpoint of the spectral evolution, only 
a small fraction of the injected energy from the Crab Pulsar needs to go to the 
magnetic field, which is consistent with previous studies. 

Subject headings: radiation mechanisms: non-thermal — ISM: supernova remnants — 
pulsars: general — ISM: individual objects(Crab Nebula) 
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INTRODUCTION 



A pulsar releases its rotational energy as a relativistic magnetized outflow called a 
pulsar wind. The pulsar wind collides w ith the surrounding super nova ejecta, forms the 



termination shock, and creates a PWN ( iKennel fc Coroniti 



1984a( ). The acceleration of 



the pulsar wind particles occurs at the terminati on shock and the PWN consists of the 



1974 



Kennel fc Coroniti 



magne tic fleld and the ultrarelativistic particles ( iRees fc GunnI 
1984bl ). Created PWN emits photons ranging from radio to TeV 7-rays via the synchrotron 
radiation and the inverse Compton scattering. Curren t status of the the o retica l models as 



well as the observational confrontations is reviewed by 



Gaensler fc Slanel toom 



The Crab Nebula is one of the best studied PWN at almost all observable wavelengths 
including its central pulsar, called the Crab Pulsar. Many studies ha ve been made to 



explain the observed properties of the Crab Nebula as a typical PWN. 



Kennel Sz Coroniti 



( 1l984a[ ) studied the spatial structure of the Crab Nebula, assuming that it is a steady 
state object (KC model). They found that the magnetization parameter a, the ratio of the 
electromagnetic energy flux to the particle energy flux just upstream the termination shock, 
must be as small as . 003 to explain the observed dynamical properties of the Crab Nebula. 



Atoyan &: Aharoniaru ( 1l996l ) succeeded to reconstruct the current observed broadband 



spectrum of the Crab Nebula by the use of the KC model. 

The KC model is not flt to consider the evolution because it is a steady state model. 
However, it is important to consider the spectral evolution of the Crab Nebula. To 
explain the flux decrea s e rates of the C r ab Nebula in radio and optical wavelengths (e.g. 



AUer fc Reynolds 



1985 



Vinyaikin 



2007 



Smith 



20031 ). we need to consider the spectral 



evolution. Moreover, we need to understand the spectral evolution of PWNe in general. 
Recent observations have found many PWNe which have a variety of characteristics in 
terms of age, expansion velocity, morphology, radiation spectrum and others. Some of these 
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characteristics may be understood by the spectral evolution of PWNe. 



For ex ample, several old PWNe which are faint in X-rays have been detected in TeV 



7-rays and Ide Jager fc Djannati-Ataj (j2008l ) discussed the possibility that some of the 
TeV 7-ray sources without an X-ray counterpart may be old PWNe. When we study the 
broadband spectrum of old PWNe, we need to take into account the evolution of the 
injected energy from the pulsar, such as the magnetic field and the relativistic particles. 
For old PWNe, most of the energy inside the PWNe is the energy which was injected when 
they were young because the spin-down power of the pulsar decreases fast with time. 



Several spectra 



de Jager et al. 



2008 



evolution moc 


els 0:' 


Gelfand et al. 


2009) 



PWNe have been studied (e.g. IZhang et al.ll2008l : 



energy injection from the pulsar, several issues still remain to be further clarified. These 
studies assume dif ferent evolution model of the particle distribution inside the PWNe. 



Zhang et al. 



(I2OO8I ) considered the escape of the particles from the PWN which has a fixed 
volume. Because a PWN is expanding while it is confined in the expanding supernova 
ejecta, it is more realistic to regar d that partic l es rat her s uffer from an adiabat ic loss 



than they escape from the nebula. 



Zhang et al. 



fl2008h and 



de Jager et al. 



(12008h used the 



conversion efficiency of the spin-down power to the particle energy ~ 10 — 50% to normalize 
the observed flux. In their model, the magnetic fleld inside the PWN does not relate to 
the spin-do wn power and the renaaining 50 — 90% of the spin-down power is not explicitly 



considered. 



Gelfand et al 



Atovan &: AharonianI (119961 ) has a similar problem of the flux normalization. 



(I2OO9I ) did not compare their calculated spectrum with the observations 



because they were interested mainly in the dynamical evolution of the PWN. 

In this paper, we revisit a spectral evolution model of PWNe, paying attention to the 
issues mentioned above. We do not take into account effects of the spatial structure of 
the PWN because it is somewhat costly and too detailed to discuss the whole radiation 
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spectrum. In our model, the PWN is simply treated as an expanding uniform sphere. The 
energy inside the PWN is injected from the pulsar spin-down energy, which is divided 
between the magnetic field and the relativistic particles with a constant ratio. Our simple 
model can describe the observed basic features. Of course, many details are not treated 
in this simplified model, such as the filamentary structures and the spatial variations of 
photon indices. 

We study the spectral evolution of the Crab Nebula as the first application of our 
model. The Crab Nebula can be used as a calibrator of our model when it will be applied 
to other PWNe in future. In section [21 we describe our model of the PWN evolution. In 
section |3l we apply this model to the Crab Nebula. In section HI discussions and conclusions 
are made. 



2. THE MODEL 

For the calculation of the spectral evolution, we need to specify the evolution of the 
magnetic field and the particle distribution function. Here, we describe the assumptions of 
our model. 



2.1. Model of the Energy Injection 



In this paper, we assume that the PWN is a uniform sphere expanding at a constant 
velocity wpwN- The assumption of the constant velocity is the easiest way to take into 
account the ex pansion of the PWN, although the real behavior of the expansion must be 



more complex. 



Gelfand et al. 



(|2009[ ) investigated the dynamical evolution of the PWN 
surrounded by supernova ejecta. They showed that the PWN expands at an almost constant 
velocity at first and then it turns to shrink and bounces in a late phase of the evolution 
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after about lOkyr from birth. We consider that the age of the PWN is younger than lOkyr 
in this paper and that the constant velocity is a good assumption in this range. The radius 
of the PWN RpwN {t) at a time t is given by 

-RpWn(^) = '*^PWN^- (1) 

For the components inside the PWN, we assume that the PWN is composed of the 
magnetic field and the relativistic electron-positron plasma, both of which are injected from 
the pulsar inside the PWN. The evolution of the spin-down power L{t) is given by 

_ n + l 

L{t) = Lo(^l + ^^ (2) 

where Lq is the initial spin-down power and Tq is the initial spin-down time. Both 
parameters are fixed, if the current pulsar period P, its time derivative P, braking index n 
and age of the pulsar t^ge are known, assuming that the moment of inertia of the pulsar is 
lO^^g -cm^. 

We divide the energy injection from the pulsar into the magnetic field energy Eb and 
the relativistic particle energy using the time independent parameter 77 (0 < 77 < 1). The 
fraction parameter 77 is the ratio of the magnetic field energy injection to the spin-down 
power, i.e., 

L{t) = E,{t) + EB{t) 

= {l-rj)-L{t)+rj-L{t). (3) 

The fraction parameter 77 in our model is similar to the magnetization parameter a in the 
KC model, although they are not the same. The magnetization parameter a is the ratio 
E-Q/Eg at the pulsar wind region immediately upstream the termination shock. On the 
other hand, the fraction parameter 77 pertains to the ratio Eb/{Eb + E^) into the PWN 
region. 
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For the particle injection, we also need to determine the injection spectrum of the 
relativistic particles. Following 



Venter fc de Jager 



( 120061 ). we assume that the injection 



spectrum of the relativistic particles Qm}{l,t) obeys a broken power-law 
^ , I Qo(t)(7/7b)-^^ for7min<7<7b, 

Qmh,t)={ (4) 
[ Qo(t)(7/7b)-^^ for 7b < 7 < 7max , 

where 7 is the Lorentz factor of the relativistic electrons and positrons. We introduce 
time independent parameters 7min, 7b, 7max, Pi and p2 which are the minimum, break and 
maximum Lorentz factors and the power-law indices at the low and high energy ranges of 
the injection spectra, respectively. Below, we use terms the low/high energy particles to 
refer to the particles which have the energy lower /higher than 7b. Note that most of PWNe 
have value pi < 2 and p2 > 2, which means that the energy of the particle injection is 
dominated by the particles 7 = 7b- 

We require that the normalization Qo{t) satisfies the following equation 

/'7max 

{l-r])L{t)= Q,^^{^,t)^m,c''d^, (5) 



mm 



where rrie and c are the mass of an electron (or positron) and the speed of light, respectively. 
From equations (111), (HI), and ([5]), the normalization parameter Qo{t) is given by 

Lo-jl-v) ( t\-^ / 7g(pi-P2) , 7r-7gra^- tTt^V' 



2.2. Model of the Evolution of the Magnetic Field 



Because the magnetic field lines are stretching and winding, it is difficult to model 
the evolution of the magnetic field B{t) in the context of the uniform PWN. We have to 
solve the relativistic magneto hydrodynamics (MHD) equations to determine the realistic 



magnetic field evolution (e.g. 



Del Zanna et al. 



2004 ). In this paper, for simplicity, we 
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assume that the magnetic field evolution is determined in the form of the magnetic field 
energy conservation, i.e., 



7]L{t')dt' 
r]Espin(t), 



(7) 



where Espmit) is the integrated spin-down energy at a time t. From equations ([T]), ([2]) and 
(I7j), the magnetic field at a time t is given by 



Bit) 



3{n - l)r]LoTo 



IPWNJ 



1 + 1 

To 



2 

n-l 



The magnetic field approximately evolves as B{t) oc t~^-^ for t > tq. Note that although 
this magnetic field evolution model may be ad hoc, its behavior is very similar to those 



(8) 



(Rees & Gunn 


1974; 


Kennel & Coroniti 


1984a; 


de Jaeer et al. 


2009; 


Venter & de Jaeer 


2006 


)• 



Rees fc Gumj (jl974f ) considered that the magnetic field evolution is determined by the 
number of turns of the central pulsar. They considered that the magnetic field in the PWN 
was built up by the winding of field lines because of the pulsar spin. The total number of 
the turns of the magnetic field line is given by 



Nit) 



2n 



-dt'. 



(9) 



This number is proportional to the magnetic fiux, which means 

OciV(t)(i?pwN(t))-'. 



(10) 



For t > To, we obtain B{t) oc so that B{t) oc t"^-^ for n = 3 and B{t) oc t'^/^ for 

n = 2.5. Our choice of the magnetic field evolution is not so different. 
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Kennel &: Coronit il ( ll984a( ) considered the steady state spatial structure of PWN. 
In the KG model, the magnetic field increases with the distance from the pulsar till the 
magnetic pressure dominates over the particle pressure. Although the assumption of the 
steady state PWN makes it difficult to compare with our model, it seems to be natural to 
regard that the magnetic energy is a constant fraction of the total energy (~ the particle 
energy), which means 

(i?(t))^(i?PWN(t))^ OC (11) 

where we include the adiabatic cooling of the total energy, so that B(t) oc for t > tq. 
Equation fill I) is alternatively interpreted as the magnetic fiux conservation. The magnetic 
field decreases more rapidly than our model, but it is still close to our model. 



In Ide Jager et al.l (120091 ). they mentioned that the spatially averaged magnetic field 
strength of the PWN decreases as B{t) oc t~^-^ in their calculation of non-relativistic MHD 
equations. This is again close to our model. 



Finally, 



Venter &: de Jagerl ( 20061) assum e d tha t the magnetic field of the PWN evolved 



Zhang et al. 



(120081 ) used this model with k = 0.5, i.e.. 



as B{t) = Bq/{1 + it/Tf) and 
B{t) oc t^^'^ for t > T = 500yr. This is somewhat different from other models presented 
above and from ours. Based on these comparisons, we regard that our model is close to the 
reality, although other choices are also possible. 



2.3. Model of the Evolution of the Particle Distribution 

We assume that the distribution of the particles in the PWN is isotropic, so that the 
particle distribution function can be easily volume integrated to be described by the energy 
distribution N{'y,t). The evolution of the particle distribution N{'y,t) is given by the 
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continuity equation in the energy space, 



-iV(7,t) + — (7(7,t)iV(7,t)) = ginj(7,t). 



We consider the coohng effects of the relativistic particles 7(7, t) including the 
synchrotron radiation, the inverse Compton scattering off the Cosmic Microwave 
Background Radiation (CMB) and the adiabatic expansion, i.e.. 



(12) 



7(7, t) = 7syn(7, t) + 7IC(7) + 7ad(7, t). 



(13) 



The inverse Compton cooling 710(7) does not depend on time because we consider that the 
target photon field is only the CMB. As the cooling effect, we do not include the inverse 
Compton scattering off the synchrotron radiation field because it never be a more important 
cooling process than the synchrotron cooling and because it costs too much computer 
power. Note that we include it in the calculation of the radiation spectrum, as described in 
section \2a 



The synchrotron cooling 7syn(7,^) is given by (e.g. iRybicki fc Lightmanl 119791 ) 

4 a, 



7syn(75 



UB{t)Y 



3 nieC 

where ctt is the Thomson cross section and ?7b(^) is the magnetic field energy density 
calculated from equation ([8]). 



(14) 



The inverse Compton cooling 710(7) is given by (e.g. 



Blumenthal fc Gould 



mm 



7ic(7) = -^-^4 r r ''^^"^'''"^ • /(g, T,).eil-q).9iq-l/A^')du,^,, (15) 

4 meC Y Jo Jo ^ini 

where h is the Planck's constant, u^^i and Uf^^ are the frequency of the CMB photons 
and that of scattered photons, respectively, ncuBi^im) is the distribution of the CMB 
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described in equation f lT9|) below, = 47/ii/ini/(meC^), q = hu^^/ {Vg{^meC^ — hufin)), 
f{q,T,) = 2glng+ (1 + 2g)(l - g) + 0.5(1 - g)(r,g)V(l + T.q), and 6 is the step function. 

Finally, the adiabatic cooling 7ad(7,^) is given by 

. / ,x "ypwN 
7ad(7,^) = —B 7TT-7 

-KpWN(,rj 

= -f (16) 

where we use equation ([T]). Note that the adiabatic cooling is independent of the expansion 
velocity t'pwN- 

As was stated in section [H we consider that it is more reasonable to treat the adiabatic 
loss rather than the escape loss. 



Zhang et al 



(|2008[ ) considered an escape of the particles 



instead of an adiabatic loss. They treated the escape of the particles based on the Bohm 
diffusion. Because higher energy particles have larger Lamor radii, the high energy particles 
suffer from the escape loss in their models. On the other hand, the adiabatic loss is the 
dominant cooling process for the low energy particles in our model, since the high energy 
particles suffer from stronger radiative loss. 



2.4. Calculation of Spectrum of the PWN 

When the evolutions of the magnetic field and the particle distribution are determined, 
we can calculate the spectral evolution of the PWN. To compare the calculated luminosity 
with the observed flux, we assume that the radiation is isotropic. The radiation processes 
which we consider in our model are the synchrotron radiation and the inverse Compton 
scattering off the CMB (IC/CMB) and the synchrotron radiation (SSC). 



The synchrotron radiation luminosity i^v,syn(^) is given by (e.g. 



Rybicki fc Lightman 



- 12 - 



19791) 



■'iy,syn 



it) 



N{j,t)P{u,^,B{t))d^, 



(17) 



where P{u,'j, B{t)) = \/3e^B(t)/{meC^)F{u/i'c{t)) is the emissivity of the synchrotron 
radiation per particle, z/c(t) = 3e-B(t)7^/(47rmcc), F{x) = x K^i2,{ii)dy with x = v/uc 
and -ft's/s being the modified Bessel function of order 5/3. 



The inverse Compton scattering luminosity Ly^^^icit) is given by (e.g. 



Blumenthal fc Gould 



1970) 



3 CTrp/lZ/fin 



iV(7,t) 



where nph(i/ini) is the distribution of the target photon fields including the CMB and the 
synchrotron radiation. The distribution of the CMB ^cmbI'^) given by the black-body 
distribution 

(19) 



'^ph(t'ini) 



■ /(g,re)ciz/ini, 



'1^ 



57r 



exp [hu/kBTcMB) - 1' 



where /cb is the Boltzmann's constant and Tcmb is the temperature of the CMB. The 
distribution of the synchrotron radiation risyniy^t) is given by 



L 



v,syn 



it) 



4:TiRl^^{t)c hv 



' u. 



(20) 



where U ~ 2.24 is the average of U{x) = (3/2) J^d/fx) In ((x + y)/{x — y))dy in a spherical 
volume and x is the ratio of the distance from the center of the PWN to the PWN radius 
Rpw^it) (c.f. lAtoyan fc Aharonianlll996l ). 



In equation (120!) . 



Zhang et al. 



(120081 ) did not use -Rpwn(^), but a new parameter rgyn 
(< -Rpwn(^)), considering that the spread of the synchrotron radiation field is smaller than 
the radius of the PWN. Because only the synchrotron photons whose frequency below 
optical wavelengths contribute to the SSC fiux and because the observed size of the optical 
synchrotron nebula is comparable to the size of the PWN, we use i?pwN(^) in equation fl20|) . 
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3. APPLICATION TO THE CRAB NEBULA 

In this section, we apply our model to the Crab Nebula as the standard calibrator of 
our model. The Crab Nebula is one of the best observed PWN in almost all observable 
wavelengths. The results of the application to the Crab Nebula will be a landmark for 
applications of our model to other PWNe. 

First, we determine the parameters to reproduce the current Crab Nebula. We can 
obtain the information about the magnet ic field and the particle distr ibution and compare 



with those of the steady state solution of 



Atoyan fc AharonianI (jl996[ ). Secondly, we show 



the spectral evolution of the Crab Nebula in our model with the use of the determined 
parameters. Through the evolutions of the magnetic field and the particle distribution, we 
explain how the behavior of the spectral evolution can be understood. We also show that 
our spectral evolution model is in a reasonably good agreement with the observations of the 
radio/optical flux decreases. Thirdly, we give a simple argument to see the dependences on 
the fraction parameter rj of the spectral evolution. Based on the general properties of PWNe 
in this simple argument, we can apply this to other PWNe in future works. Finally, we 
discuss about the fitted parameters other than the fraction parameter rj, which characterize 
the injection spectrum (equation (j4])). 

The Crab pulsar has the period 3.31 x lO^^s, its time derivative 4.21 x 10^^'^s ■ s^^ and 
braking index 2.51. The progenitor supernova is SN1054, which means the age of the Crab 
Nebula t^ge ~ 950yr. Because all the necessary pulsar parameters are known, the evolution 
of the spin-down power in equation ([2]) is fixed, as tq ~ 700yr and Lq = 3.4 x lO^^ergs ■ s^^ 
The distance to the Crab Nebula 2kpc is used to convert the luminosity into the isotropic 
flux. The Crab synchrotron Nebula is roughly an elliptical shape with a major axis of 
4.4pc and a minor axis of 2.9pc. Here, we regard that the Crab Nebula is a sphere of the 
diameter ~ 3.5pc. Combining with tage ~ 950yr, the constant expansion velocity becomes 
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fpwN ~ 1800km/s, which is close to the observed expansion velocity of the Crab Nebula. 



3.1. Current Spectrum of the Crab Nebula 

Figure [1] shows calculated current spectrum of the Crab Nebula with the current 
observational data. The adopted parameters are shown in Table [TJ As seen in Figure [H the 
SSC flux is stronger than the IC/CMB flux in 7-rays. Note that the IC/CMB flux is almost 
in the Thomson regime, while the SSC flux is largely affected by the Klein-Nishina effect. 

We fit the data with the parameters rj = 0.005, 7max = 7.0 x 10^, 7b = 6.0 x 10^, 
7mm = 1-0 X 10^, pi = 1.5, and p2 = 2.5. The fraction parameter i] governs the absolute 
values of the fluxes and the flux ratio of the inverse Compton scattering to the synchrotron 
radiation, as discussed in more detail in section 13.31 The KC model derived a ^ 1 from 
the viewpoint of the current dynamical structure of the Crab Nebula, while we determine 
1] 1 from the viewpoint of the spectral evolution. The parameters 7max5 Ih, Pi, and p2 are 
fixed to reproduce the observed synchrotron spectral shape, such as the spectral breaks and 
the photon indices, while 7min should be regarded as an upper limit to reproduce the radio 



flux at the lowest frequency. In section 13.41 these fitted parameters characterizing particle 
injection are discussed in detail. 

In our calculation, the current magnetic field streng th of the Crab Nebula turns out to 



be -Bnow = 85/iG, which is smaller than ~ 300yuG used by 



Atoyan &: Aharonian 



differe nce of the magnetic field strength can be explained as follows. 



(1996). This 



Atoyan fc Aharonian 



(1l996l ) adopted Skc ~ 300/iG from the KC model and adjusted the particle number to 
reproduce the observations. They applied roughly half a spin-down power compared with 
the KC model to reproduce the spectrum and thus the other half is missing. On the other 
hand, all the injected spin-down power is divided between the magnetic field and the 
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particle energies in our model. If we adopt -Bnow = -Bkc ~ 300yuG, the synchrotron flux and 
also the SSC flux increase by ab out an order of magnitude. Note that the relativistic MHD 



simulation by 



Volpi et al. 



( 120081 ) also indicates a smaller value of the spatially averaged 
magnetic field strength ~ lOO/iG, which is close to our value -Bnow = 85/iG. 



3.2. Spectral Evolution of the Crab Nebula 

Our model can calculate the past and the future spectra of the Crab Nebula. All 
the calculated results in this section shown in Figures [2ll5] are with the use of the same 
parameters as Figure [TJ 

Figure [2] shows the spectral evolution. As seen in Figure [21 the synchrotron flux 
decreases with time and the SSC flux also decreases with time in accordance with the 
synchrotron flux, while the IC/CMB flux decreases more slowly than the SSC flux. We 
discuss their time dependence in section 13.31 An important result in Figure [2] is that the 
flux ratio of 7-rays to X-rays increases with time. This supports the view that old PWNe 
can be observed as 7-ray sources with no or a weak X-ray counterpart. 

The radio/optical observations of the Crab Nebula have suggested that the radio/optical 
flux of the Crab Nebula is decreasing with time. The inferred rate of the radio flux decrease 
i s — 0.17 ± 0.02% /yr and almost independent of the frequency for a range 86 - 8000MHz 



( IVinyaikin 



20071 ). Our model predicts the current rate ~ — 0.16%/yr, which is almost 



consistent with the observation. T he inferred ra te of the optical continuum flux decrease 



is -0.55%/yr calibrated at 5000A ([Smith 



20031). Our model predicts the current rate 



~ — 0.24%/yr, which is a factor of two smaller than the observation. However the trend that 
the decreasing rate increases with frequency matches the observations. This is because the 
optical emitting particles suffer from stronger synchrotron cooling than the radio emitting 



- 16 - 



particles. 

To understand the detailed features shown in Figure [21 we examine features of the 
evolution of the particle energy distribution in Figure |3l For comparison, the injected 
spectrum of the particles till lOkyr without the cooling effects is plotted by the dot-dashed 
line in Figure El The evolution of the particle distribution is characterized differently in four 
energy ranges. First, for 7 > 10^, the particle number increases with time. This increase 
of the high energy particles relates to an increase of the IC/CMB flux in Figure [2l after 
Ikyr above lOTeV. Secondly, for 10^ < 7 < 10^, the behavior of the particle distribution 
is complex. The particle number and the power-law index in this energy range do not 
monotonically change. The radiation spectrum in the range from infrared to X-rays reflects 
this complex behavior. Thirdly, for 10^ < 7 < 10^, the change of the particle distribution 
is small. The difference between all the lines for several epochs is less than a factor of two. 
This leads to an important conclusion that the radio flux decrease is mainly due to the 
decrease of the magnetic field, so that the observations of the radio flux decrease support 
our model of the magnetic field evolution. Lastly, for 7 < 10^, there exist particles with 
7 < 7min- This is because the adiabatic cooling is still effective at low energy. We discuss 
more about the features of the particle 7 > 10^ in later paragraphs, which relate to the 
spectral evolution above lO^Hz. As discussed below, these features indicate the importance 
to consider the evolution of the particle injection rate and the cooling time due to various 
processes. 

FigureHlplots the evolution of cooling time r(7, t) = 7/17(7, t)| and helps to understand 
features of the particles with 7 > 10^ in Figure [31 As seen in Figure [H at an age of Ikyr, the 
synchrotron cooling is effective for the particles with 7 > 10^, while the adiabatic cooling is 
effective with 7 < 10®. However, because of a rapid decrease of the magnetic field energy 
density, the adiabatic cooling dominates at the 7 < 10* at an age of lOkyr. From equations 
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dH]), dHl) and f|T6l) . Tgyn oc 7 for t < tq and Tgyn oc 7 ^t^ for t > tq, while Tad oc t at any 
time. 

We examine the increase of the particle number for 7 > 10^ shown in Figure 
|3l The increase means that the particle injection dominates the cooling effect in this 
energy range. Equation (IT^ is approxmately expressed as {d/dt ~ t^^ = and 

iV(7,t)~^^^^^^-Qinj(7,t)- (21) 

"^01711 ~r Tcool 

Because synchrotron cooling dominates in this energy range Tcooi ~ Tsyn < Tdyn, equation 
fl2T|) becomes ~ TgynQinj- For t < tq, the injection term behaves as Qinj oc 7~^^, so that 
oc 7~^2~^t^. For t > tq, the injection term behaves as Qinj oc ^~P2j(-7/3 ^ _ 2 5^ gg 
that N (X ^y' 

~P2-i^2/3^ This is the reason why the particle number increases with time and 
the particle distribution is softer than the injection distribution for 7 > 10*. 

We next examine the steadiness of the particles number for 10^ < 7 < 10^ shown in 
Figure [3l For t < tq, because TcooI ~ fad ~ Tdyn in equation (1211) . we have oc •y'^^t. For 
t > Tq, because the injection term more rapidly decreases with time than the adiabatic 
cooling term, equation (fT2|) is approximated by 

^iV(7,t)~|-f>(7,t)). (22) 



at ' (97 Vt 

When we assume a form oc 7~"t^, we obtain P = 1 — a. For a = pi = 1.5, we have 
P = —0.5. The particle number increases until t ~ tq, and then it decreases slowly with 
time as A^ oc t'^-^ for 10^ < 7 < 10^ 

The feature of the particle distribution for 10^ < 7 < 10* shown in Figure [3] is complex 
because the dominant cooling process changes with time as seen in Figure HI Until an 
age of a few thousand years, the feature is similar to that for 7 > 10* (equation fl2T]) with 
Tcooi ~ Tsyn), i-G., the partlclc number increases with time and the particle distribution is 
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softer than the injection distribution. After that, the feature becomes similar to that for 
10^ < 7 < 10^ for t > tq (equation f l2^ ) and the particle number decreases with time. 
Because the time when r^d ~ Tgyn is different for different energy, the particle distribution 
becomes harder than oc 7~P2-i at an age of lOkyr. 

In Figure we show the evolution of each energy component: the particle energy, the 
magnetic field energy, the radiated energies via the synchrotron radiation and the inverse 
Compton scattering and the wasted energy via the adiabatic expansion, all of which are 
normalized by the integrated spin-down energy Espin(t). As seen in Figure [5], although 
almost all the injected energy goes to the particle energy, the particle energy decreases with 
time by the cooling effects. The magnetic field energy is always constant fraction of the 
injected energy rjEspmit) as assumed in equation ([7]). Most of the particle energy is cooled 
by the synchrotron radiation in the early phase (< Ikyr) because the cooling break due to 
the synchrotron cooling is a little smaller than 7b as seen in Figure [31 On the other hand, 
most of the injected energy goes to the adiabatic loss in the late phase (> Ikyr), which 
means that most of the injected energy goes to the kinetic energy of the supernova ejecta. 
The inverse Compton cooling is not important in the case of the Crab Nebula for t < lOkyr, 
which is also seen in Figure |H Note that the integrated spin-down energy Espin{t) also 
increases with time, for example, i?spin(100yr) = 9.1 x lO'^^ergs, i?spin(lkyr) = 3.9 x lO^^ergs 
and -Espin(lOkyr) = 5.5 x 10"^^ ergs. 



3.3. Characteristics of the Spectral Evolution 

We discuss the evolution of powers of the synchrotron radiation and the inverse 
Compton scattering and their dependence on the fraction parameter rj, which will help to 
apply our model to other PWNe in future. 
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We assume t > tq in this section aiming to discuss primarily old PWNe. For simplicity, 
we also assume that the particle energy is expressed as N{'y^,t)'y'^meC^ = (1 — ri)C,it)Espi^it), 
where ^{t) < 1 accounts for the cooling effects. As seen in Figure A^(7b, t)7b^eC^ ~ const, 
and ^(t) ~ const. ~ 0.1 is a fairly good approximation. Moreover, we assume that powers 
of all radiation mechanisms are dominated by the emission from the particles 7 = 7b- 
This assumption is valid when the low energy power-law index < 2 and the high energy 
power-law index > 3 in the particle distribution. 

The power ratio of the inverse Compton scattering to the synchrotron radiation is 
given by 



Picit) , f/ph(t) 

~ /KN- 



(23) 



Psynit) U^{t) ' 

where /kn < 1 represents the Klein- Nishina effect and f/ph(t) is the energy density of the 
target photon field, /kn ~ 1 for the IC/CMB and /kn < 0.1 for the SSC. 



The power of the synchrotron radiation is given by (e.g. 



Rybicki Sz Lightman 



19791) 



Psjnit) ~ ^aTC7^f/B(t)-7biV(7b,t) 

where 7bA^(7b,^) is the number of the particles around 7 = 7b and we use equation 
([7]). Note that although the power of the synchrotron radiation has been conventionally 
compared with the spin-down power L{t), the power of the synchrotron radiation relates to 
the integrated spin-down energy Espin{t) rather than the instantaneous spin-down power 
L(t) in equation (!24l) . 

Now equation fl23|) becomes, in the case of the IC/CMB (f/ph(i) = Uqub), 
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Pic/CMB{t) „ /^^WVrr 
lUtT ^ ^™l^J """"" 

oc %^ocrV. (25) 
In the case of the SSC, we estimate the synchrotron radiation energy density Usyn{t) as 



Then we obtain 



^g^T--^^- / ,-^-^oc(l-r,)t-^ (27) 

where rsyn(7b,t) = A^(7b, ^)7b"^eC^/-Psyn(i) is the synchrotron coohng time of the particles 
with 7 = 7b at time t. 

All the powers Psyn(t) oc r/(l - r])t-^, Pic/cMB(t) oc (1 - r/), and Pssc(^) oc r/(l - r])H-^ 
depend on the fraction parameter t] and time in different ways. These general characteristics 
of the radiation powers make it possible to estimate the fraction parameter r] from the 
viewpoint of the spectral evolution. For example, the time dependence of Pssc(^) and 
Pic/CMB(i^) can explain that the SSC flux is dominant in high energy 7-rays in Figure [H 
Using the same parameters as in Figure [H we find Pic/CMB(^age)/Psyn(^age) ~ 10^'^ and 

Pssc(^age)/Psyn(^age) ~ 10~^. 



In 



de Jager fc Djannati-Atai 



(120081 ) , they suggested that some of the Te V 7-ray sources 
without an X-ray counterpart can be old PWNe because the X-ray emitting particles are 
cooled more rapidly than the TeV 7-ray emitting particles by the synchrotron cooling. In 
our model, old PWNe can also be the 7-ray sources with no or weak X-ray counterpart 
because the IC/CMB is almost time independent. However the same result comes from the 
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different reason. Because the number of the high energy particles is even increasing in our 
model as seen in Figure El this result comes from the rapid decrease of the magnetic field. 
PWNe which are old and has a small fraction parameter rj can be recognized as the 7-ray 
source without an X-ray counterpart. Note that -Pic/cmb(^) slowly decreases with time, 
when we properly consider the cooling effects C,{t), as seen in Figure [2] at 3kyr and lOkyr. 



3.4. Parameters Characterizing Particle Injection 

The fitted parameters pi, p2, 7max, 7b, and 7min relate to the broken power-law injection 
in equation (|4]). The above parameters include the information about the acceleration at 
the termination shock and the physics of the pulsar wind and the pulsar magnetosphere. 
We discuss about them in the framework of our model. 

For the low and high energy power-law indices at injection pi = 1.5 and p2 = 2.5, we 
determine these values to reproduce the observed radio and X-ray photon indices, Fr and 
Fx. Radio photon index F^ is related to pi, but X-ray photon index Fx is not simply related 
to P2 because of the synchrotron cooling. Because it is difficult to make a large spe ctral 



Reynolds 



20091), we 



break Fx — F^ = AF > 0.5 from the single power-law injection (e.g. 
adopt the broken power-law injection. 

For the maximum energy 7max = 7.0 x 10^, we determine this value to reproduce the 
observed spectral break at ~ lOOMeV. The maximum energy 7max is conventionally given 
by comparing the acceleration time with the cooling time or the size of the acceleration site 
with the Lamor radius. Both conditions give 7max ~ 10^'' using the magnetic field strength 
~ 100/iG and the size of the acceleration site ~ O.lpc. The fitted value is near the limit of 
theoretical expectation. 

For the break energy 7b = 6.0 x 10^, we determine this value to reproduce the observed 
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spectral break around optical wavelengths. Although the KC model related 7b to the bulk 
Lorentz factor of the pulsar wind immediately upstream the termination shock, our model 
does not allow this connection as discussed below. 

Finally, for the minimum energy 7min = 1.0 x 10^, we determine this value to 
reproduce the flux of the observed minimum frequency at radio wavelengths. Because 
Pi > 1, the particles around 7min determine the injection of the particle number as 
Ne{t) ~ <5inj(7min, ^)7min- The particle number conservation outside the pulsar light 
cylinder leads to the particle number flux into the PWN is much larger than the 
Goldreich- Julian number flux ncj, which is the particle number flux from the pulsar 
polar cap. In our model, the pair production multiphcity at an age of a thousand year 
K ~ (-^(t)/7b"^eC^)(7b/7min)^^~^ ~ 10^ from equations (jlj) and ([5]). Theoretically, the 
multiplicity is estimated as k ~ 10^ — 10^. Our value n ~ 10^ is somewhat large. The mean 
energy of the injected particles L(t)/Ne(t) ~ 'Jh^ec'^ilh/lmm)'^^'^^ is significantly smaller 
than 7b meC^. 



Atoyan fc AharonianI (Il996l ) did not use the broken power-law injection. They divided 
the particles inside the PWN into the low energy particles as the relic electrons and the high 
energy particles as the wind electrons. This assumption can reduce the multiplicity, but 
the origin of the relic electrons becomes another problem. It is difficult to discuss anything 
more about this problem from our spectral evolution model. One thing what we should 
note is that the relation between the radio flux decrease and the magnetic field evolution is 
always kept, i.e., the relation is independent of how and when the radio emitting particles 
are injected. This is because the low energy particles whose power-law index is pi do 
exist from the observation and their distribution hardly changes by the cooling effect, as 
discussed in section [321 
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4. DISCUSSIONS AND CONCLUSIONS 

The evolutions of the magnetic field and the particle distribution determine the 
spectral evolution of a PWN. The evolution of the particle distribution is affected by 
the assumptions of the magnetic field evolution model, uniformity of the PWN, particle 
injection spectrum, and expansion evolution of the PWN. We discuss about the effects of 
these assumptions which are made in our model. Finally, we summarize the conclusions of 
this paper. 

4.1. Discussion 

Our model of the magnetic field evolution is somewhat ad hoc. As discussed in section 
12.21 however, the time dependence of the magnetic field strength B oc t^"' for t > tq seems 
to be in a range of 1.5 < a < 2.0 from other theoretical considerations and we adopt 
a = 1.5. Moreover, because our result of the radio fiux decrease of the Crab Nebula is 
almost consistent with the observation, our model of the magnetic field evolution can be 
near the truth. 

For the assumption of the uniform PWN, many non-uniform structures have been 
observed, such as the filamentary structures and the spatial variations of photon indices. 
However, for the calculation of the total spectrum of the PWN, the energetics of the PWN 
is important in the lowest order. We consider that the assumption of the uniform PWN is 
reasonable for the calculation of the total spectral evolution of the PWN. 

For the injection spectrum of the particle distribution, the acceleration of the particles 
is an unsolved problem and we adopt the broken power-law injection. It should be noted 
that one of the important conclusions in our study that old PWNe can be observed as 
7-ray sources with no or weak X-ray counterpart is not affected by the broken power-law 
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assumption. This is because low energy particles do not contribute to X-ray and high 
energy 7-ray emissions. 

The use of the time independent parameters 7min, 7b and 7max can be improved as 
time dependent parameters, because the physical condition of pulsar wind termination 
shock may be change with the decr ease of the spin-down pow er of the pulsar. Considering 



the time dependence of 7max(^), in 



Venter fc de Jagerl (120061 ) . they use the condition 



ri^{t) < O.Srg with time dependent magnetic field, where tl is the electrons' Larmor radius 
and Ts is the radius of the termination shock. As discussed in section 13.41 our model 
satisfy this condition. Both 7jnin and 7b are important parameters because these may 
include the information about the pulsar magnetosphere and the pulsar wind. However, 
the time dependences of them are uncertain. For simplicity, we used all of them as the time 
independent parameters in the present paper. 

Constant velocity expansion is a good assumption for young PWNe, although the 



expansion of t he PWN shou 



the PWN (e.g. 



Gelfand et al. 



d be calculated by taking account for the environment of 



20091 ). In our model, the magnetic field evolution explicitly 
depends on the expansion velocity fpwN (see equation ([8])). To understand a little more 
about the effects of the expansion evolution, we study how the Crab Nebula would be 
observed in the context of the constant velocity expansion. The Crab Nebula is one of 
the sources without observable SNR shell. It may be because the surrounding interstellar 
medium is less dense than other PWNe with observable SNR shell. If the Crab Nebula were 
in the different surroundings, the expansion velocity would also be affected. That is, if it 
were in a less or more dense surroundings, the expansion velocity would be more rapid or 
slower. 

An example of a twice rapid expansion is shown in Figure El All the parameters except 
the expansion velocity are the same as in Figure [H An example of a half velocity expansion 
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is shown in Figure [7] and all parameters except for the expansion velocity are the same 

as in Figure [TJ Note that both of them are hypothetical PWNe, not the Crab Nebula 

itself. For the spectrum of the rapid expansion case, the absolute value of the flux becomes 

smaller and the flux ratio of the inverse Compton scattering to the synchrotron radiation 

is larger than the real Crab Nebula shown in Figure [T] and for the spectrum of the slow 

expansion case vice versa. Comparing the particle distribution in Figure with that in 

Figure Ul the low energy particles take the same distribution, but the high energy particle 

distribution in Figure [7] is steeper than that in Figure [61 These spectral behaviors against 

the expansion velocity can be understood from equations ([H]), ([HD and f[T^ . In our model, 

the magnetic field becomes smaller when the expansion velocity becomes larger. This leads 

to the difference in the absolute flux and the synchrotron cooling which changes the high 

energy particle distribution. On the other hand, because the adiabatic cooling does not 

depend on the absolute value of the expansion velocity, the low energy particle distribution 

does not change. 

I 



Lastly, in lAtoyan fc AharonianI (119961 ). they included the infrared photons and the 



starlight for the target photon fields of the inverse Compton scattering. Alth ough these 



soft p hotons can significantly contribute to the 7-ray flux of other PWNe (e.g. 



Porter et al. 



20061 ). this is not the case of the Crab Nebula. Because the Crab Nebula is located far 



away from the galactic center (~ lOkpc) and galactic plane (~ 200pc), the energy density 
of these soft photon fields is less than the solar neighborhood (~ 8kpc from the galactic 
center). Even if we assume that the energy density of the infrared photons and the starlight 
is the same as the solar neighborhood, the inverse Compton scattering off these photon 
fields contributes less than 30 % of the current total 7-ray flux. Note that it also does not 
much affect the 7-ray spectrum, when the SSC flux decreases if the energy density of these 
soft photon fields is less than the half that in the solar neighborhood since inverse Compton 
scattering off the CMB dominates there. 
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4.2. Conclusions 



In this paper, we built a model of the spectral evolution of PWNe and applied to 
the Crab Nebula as a calibrator of our model. We solved the equation for the particle 
distribution function considering adiabatic and radiative losses with a simple model of 
magnetic field evolution. 

The flux decrease of the 7-rays is more moderate than radio to X-rays, because the 
magnetic field decreases rapidly, which implies that ol d PWNe can be observed as 7-ra y 
sources with no or weak X-ray counterpart. Although 
obtained the same result but for a different reason that the X-ray emitting particles are 
cooled more rapidly than TeV 7-ray emitting particles. 

The current observed spectrum of the Crab Nebula is reconstructed when the fraction 
parameter has a small value 77 = 0.005. Thi s is consistent with the pr ediction of the 



de Jager fc Djannati-Atail ( 120081 ) 



magnetization parameter a ^ 1 obtained by 



Kennel fc Coroniti 



f!1984al ). They obtained 



o" ^ 1 from the viewpoint of the current dynamical structure of the Crab Nebula, while we 
determine 77 ^ 1 from the viewpoint of the spectral evolution. 



The smaller value of the current magnetic field 



85/iG is nee ded to recoii s truct 



the observed spectrum of the Crab Nebula. This is consistent with that of 



Volpi et al. 



(120081 ) 



for the spatially averaged magnetic field strength ~ lOO/xG from th eir relativistic MHD 



simul ation, but smaller than ~ 300/iG in most of other papers (e.g. 
19961). 



Atoyan fc Aharonian 



Our model can predict the spectral evolution of the Crab Nebula, and the observed 
flux decrease of the Crab Nebula at radio wavelengths can be explained by our model. This 
conclusion does not depend on the assumption of the broken power-law injection, and gives 
the validity of the our magnetic field evolution model. The observed flux decrease of the 
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Crab Nebula in optical wavelengths is somewhat larger than our model, but the trend that 
the decreasing rate increases with frequency matches observations. 

The minimum energy 7inin is related to the pair production multiplicity in the pulsar 
magnetosphere, since low energy particles are assumed to be injected in the same way as 
high energy particles in our model. Our result of the minimum energy 7i„in = 1.0 x 10^ and 
the low energy power- law index pi — 1.5 means that the multiplicity k ~ 10^ is necessarily 
larger than other models which adopt a separate origin of low energy particles. 

We are grateful to Y. Ohira for useful discussions. This work is partly supported by 
KAKENHI (F. T. , 20540231) 
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Fig. 1. — Current spectrum of the Crab Nebula in our model and the observational data. 
The solid line is the total spectrum which is the sum of the synchrotron (dotted line), 
IC/CMB ( dashed line), and S SC (dot- dashed line) spectra , respe ctively. The obs erved data 



taken from 



fll979h: 



Baars et a. 



Temin et al. 



Aharonian et al. 



1977 ) (radio). 



(I2nn6h: 



(12004 . 



20061); 



Macfas-Perez et al. 



Nev fc SteinI f 



Albert et al 



9681) ( IB). 



(120081 ) 



(120101) f rad i o-opt ical) , 



Grasdalen 



Kuiper et al. 



Abdo et al. 



jpiD (X-ray-7-ray), 
20101 ) (very high energy 



7-ray). Used parameters are tabulated in Tabl^Il 
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Fig. 2. — Spectral evolution of the Crab Nebula. The thin solid line is 300yr from the birth. 
The thick solid, thin dotted and thin dashed lines are Ikyr, 3kyr, and lOkyr from birth, 
respectively. Each line represents the total spectra which are the sum of the synchrotron, 
IC/CMB and SSC spectra. Used parameters are the same as in Figured] 
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Fig. 3. — Evolution of the particle distribution. The thin solid line is the distribution at 
300yr from the birth. The thick solid, thin dotted and thin dashed lines are those at Ikyr, 
3kyr, and lOkyr, respectively. The dot-dashed line is the total injected particles at an age 
of lOkyr. Used parameters are the same as in Figure [TJ 
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Fig. 4. — Cooling times as a function of the Lorentz factor. The thin and thick sohd hnes are 
the total cooling time rcooi(7,^) at t = Ikyr and lOkyr, respectively. The dotted, dashed and 
dot-dashed lines are rsyn(7, lOkyr), rad(7, lOkyr), and ric(7), respectively. Used parameters 
are the same as in Figure [1] 
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Fig. 5. — Evolution of the energy content inside the Crab Nebula and the radiated en- 
ergy and wasted energy by adiabatic expansion. The thin solid line corresponds to the 
particle energy P'^'^'' 'ymcc'^N('y,t)d'y, the thick solid line is the magnetic field energy 

^ Tmin 

{4:n/3)RpT^j^(t)UB(t), the dotted, dashed and dot-dashed lines are the radiated energy via 
synchrotron radiation J^'"'"' /J mgC^ |7syn(7, ^OMT'^^'j the wasted energy via adiabatic expan- 
sion n™'"' f!^ meC^\'yad('l,t')\d'ydt' , and the radiated energy via inverse Compton scattering 



J^™''" mec'^\'yic{'y)\d'ydt' respectively. All the lines are normalized by the integrated spin- 



7min Jo 

down energy i?spm(i)- Used parameters are the same as in Figure [H 
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Fig. 6. — Evolution of the emission spectrum (left panel) and particle distribution (right 
panel) of the PWN for the rapid expansion case. Used parameters are the same as in Figure 
[1] except for the expansion velocity being twice (fpwN = 3600km/s). 
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Fig. 7. — Evolution of the emission spectrum (left panel) and particle distribution (right 
panel) of the PWN for the slow expansion case. Used parameters are the same as in Figure 
[1] except for the expansion velocity being a half (fpwN = 900km/s). 
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Table 1: Used parameters to reproduce the current observed spectrum of the Crab Nebula. 



A(l()])to(l Pai-anieter 


Syinl)()l 




Current Period (s) 


P 


3.31 X 10-2 


Current Period Derivative (s • s~^) 


P 


4.21 X 10"^^ 


Braking Index 


n 


2.51 


Age (yr) 


4ige 


950 


Expansion Velocity (km/s) 


fpWN 


1800 


Fitted Parameter 


Fraction Parameter {E^/{Eq + E^)) 


V 


0.005 


Low Energy Power-law Index at Injection 


Pi 


1.5 


High Energy Power-law Index at Injection 


P2 


2.5 


Maximum Energy at Injection 




7.0 X 10^ 


Break Energy at Injection 


7b 


6.0 X 10^ 


Minimum Energy at Injection 




1.0 X 10^ 



